Interaction of quantum Hall systems with waveguide elastic modes 
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An interaction of non-uniform plane elastic modes of the waveguide type with monolayer and 
double-layer quantum Hall systems is considered. It is shown, that unlike the case of the surface 
acoustic wave propagation, the restriction on maximal values of the wave vectors for which the 
velocity shift can be observed experimentally does not take place for the waveguide modes. In case 
of study of incompressible fractional quantum Hall states, the effect can be used for measuring a 
dependence of the effective magnetic length on the filling factor and for observing phase transitions 
in double-layer systems under the interlayer distance variation 
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I. INTRODUCTION 



The interaction of surface acoustic waves (SAW) with 
quantum Hall systems was intensively studied both ex- 
perimentally [|],[| and theoretically fj|| during the last 
few years. Measurements of both velocity shift and ab- 
sorption of SAW's provide an indirect method of studying 
the dynamical conductivity of a two dimensional clcc- 
tron(2DE) gas. The SAW experiments (together with a 
study of the magnetic focussing effect and of a 

temperature dependence of the d.c. conductivity of the 
quantum Hall system Jl^Tifl) are considered as an exper- 
imental proof for the composite fermion model ]l5| , ^6| ,p[ 
of the fractional quantum Hall(FQH) effect. 

In experiments, the SAW is localized on the sur- 
face of the sample, while the two dimensional electron 
layer is placed at some depth d (typical values of d are 
10 3 -i- 5 • 10 3 A). Therefore the study of a momentum 
(q) dependence of the conductivity is restricted by the 
inequality q < qo, where go is of order of d . 

The measurements beyond that region are especially 
important for an investigation of dynamical properties of 
incompressible FQH states. It was shown in Ref. [ p"7| |, 
that at filling factors (v) corresponding to such states an 
interaction of acoustic phonon with 2DEs results in an 
essential phase velocity shift (Av) for phonons at finite 
wave vectors. According to Ref. |17 , the function Av(q) 
oscillates at finite q. The period of oscillation is ~ K-ff^ 
where A e // is the effective magnetic length. Therefore, 



the measurements of Av(q) can be used as a way of an ob- 
serving the X e f f (v) dependence. It was shown Ref. Q , 
that the measurements of the velocity shifts of high fre- 
quency phonons interacting with a double layer FQH 
system can be also used as a method of a registration 
of phase transitions between FQH states described by 
different quantum numbers. (They are the generalized 
Laughlin states, which correspond to the Halperin wave 
function |l9| ] . The composite fermion theory for such 
states was developed by Lopez and Fradkin [p0|] .) 

The values of the velocity shifts obtained in 



Refs. [[L7 18 1 correspond to bulk acoustic waves interact- 
ing with a superlattice of 2DE layers. For a single mono- 
layer (or a double layer) system, the elastic wave should 
be localized near the electron layer to provide the values 
of the effect, which can be observed experimentally. At 
large q it is not the case of SAW, which is localized near 
the surface. The experimental geometry required can be 
realized in layered elastic media, where waveguide elastic 
modes may exist. We consider the case of an infinite elas- 
tic medium which incorporates an elastic layer with the 
acoustic properties different from the ones for the bulk. 
The interaction of non-uniform transverse elastic waves 
propagating in such a medium with a quantum Hall sys- 
tem is studied. We consider the case, for which the waves 
vectors and the polarization vectors of the non-uniform 
modes are parallel to the interface boundaries. The ve- 
locity of the transverse sound for the layer is assumed 
to be smaller then for the bulk matrix (it is the condi- 
tion of the existence of such modes). The piezoelectric 
interaction mechanism is implied. 

Heterostructures Al x Ga\~ x As — GaAs used in exper- 
imental studies of double layer quantum Hall systems 
obey the elastic properties required. For instance, the 
structures A - B - A- B - A and A - B - A (where A 
- Al x Gai^ x As, B - GaAs) were used in Refs. [§l|,|2j. 
Such structures can be approximately described by a 
model of a bulk sample of AlGaAs incorporating a layer 
of GaAs. (One can expect that the effect of the sur- 
face, the substrate and the thin central layer (for the 
A — B — A — B — A structure) can be neglected for the 
wave vectors considered). The transverse sound velocity 
for GaAs is smaller then for Al x Gai_ x As. Therefore, 
these exist non-uniform elastic modes localized near the 
GaAs layer. As it follows from our results, the values 
of the velocity shifts of those modes are within the ex- 
perimental resolution, while they are smaller than in the 
SAW experiments. We also consider the geometry for 
which the heterostructure is places in the matrix of an- 
other material. For the last case the values of the effect 
are of the same order or higher than for the SAW. 
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II. MATRIX ELEMENTS OF INTERACTION 

Let us start from a simple problem of renormalization 
of the phase velocity of acoustic phonons interacting with 
2DEs. The phonon part of Hamiltonian has the form 



ph 



(i) 



the expressions obtained. This choice corresponds ap- 
proximately to the parameters for AlGaAs and GaAs.) 
We consider the case when the axis x is chosen along the 
[100] direction and the z axis - along the [001] direction 
Let the wave vector and the displacement vector of the 
elastic wave are directed along the x and y axes corre- 
spondingly. The wave equation for the Fourier compo- 
nent of the displacement u y (q, z,u>) = u reads as 



where tu q = vq is the phonon spectrum, b + (b) are the 
phonon creation (annihilation) operators. The Hamilto- 
nian of the electron-phonon interaction is chosen in the 
form 



(2) 



where ^ ,+ (^E') are the electron creation (annihilation) op- 
erators, k is the electron layer index, S is the area of the 
layer, g are the matrix elements of the interaction. 

The velocity shift under such a choice of H v h Hint is 
described by the expression 



— = —y29q k 9qk'D kkl (q,vq), 
v vq f—' H 

k,k 



(3) 



where D{q,u>) is the electron density-density response 
function. In the composite fermion approach p|,p!6| an 
additional statistical interaction between electrons is in- 
troduced, which results in a modification of the functions 
D in Eq. (§|). 

We will use the approach similar to that of Ref. Q] for 
the calculations of the matrix elements g qk - The piezo- 
electric interaction of the elastic wave with the electrons 
corresponds to the Hamiltonian 



H 



E 



r etprk 



(4) 



where <p is the electric potential, generated by the elas- 
tic mode. According to the linear response theory tp in 
Eq. (U) is the external potential. Therefore, to calculate 
its value one should neglect the influence of the electron 
system. 

We will use the smallness of the piezoelectric interac- 
tion (the explicit expression for the small parameter is 
given below). Let us discuss the solution of the elastic 
problem with the piezoelectric constant set to zero. Let 
the interfaces be in the x-y planes with z = ±a coor- 
dinates. We denote q, the transverse sound velocity in 
the layer (\z\ < a), and q,, the transverse sound velocity 
in the bulk matrix (\z\ > a). The condition q, > q is 
implied. For definiteness, let us assume that the crystal 
symmetry of the bulk matrix and of the layer is the cubic 
one, and the elastic moduli C44 are the same for the layer 
and for the matrix, while they have different densities. 
(We use such a special choice of parameters to simplify 



/ 2 2 2\ n 
g^-( C l(b)'l -W )U = 0. 



(5) 



One can find from the boundary conditions that the func- 
tion u(q,z) and du(q, z)/dz are continuous ones at the 
interfaces. The localized solutions of Eq. (||) have the 
form 

u e (q, Z ) = A q \^f- XqilZl - a) \ Z \^ a (6) 
y ' q \ cos(r]qz) \z\ < a 



(the even mode), 



sin(rjqa)e~ Xq( - z ~ a) z>a 
u a (q, z) = A° ^ sm(rjqz) \z\ < a (7) 

- sin(?7<ja)e A,3 ^ +a) z < -a 



(the odd mode). In Eqs. 



I) A = VI -v 2 /ci, r, = 



\J v 2 /cf — 1. The velocity v for the even mode is given 
by the solution of the equation 

t&n(r]qa) = A/77, (8) 

and for the odd mode - by the solution of the equation 

cot(r]qa) = —A/7/. (9) 



At q < 7r/ (2ay/c^/cf — 1) there exists only one (even) 
mode. At larger q other solutions of Eqs.(||J^) emerge 
(waveguide modes). 

The transverse component of the displacement vector 
can be rewritten in terms of b operators: 



u e (o)(r,z) = ^2u e{o) (q, z)e 2qr (b q + bt q )- 



(10) 



Inserting u e i Q \ (r, z) into the expression for the energy of 
elastic vibrations 



E u = 1 J d 2 rdz[p{z) {du/dtf 
+C44 ((du/dx) 2 + (du/dz) 2 



(11) 



and comparing the result with Eq. (|l|), we determine the 
values of the normalization factors in Eqs. (MW): 



A--A°-- I V 
q ~ q ~ aV2Sc 44 ' 



(12) 



where a = (qa(l + of ) + 1/A) 1 / 2 . 

The electric potential can be found from the solution 
of the Poisson's equation 
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AyJ = -—pi,jmdiUj m , 



(13) 



where Uj m is the strain tensor 



e is the dielectric con- 
stant, Pij m is the piezoelectric tensor. In cubic crystals 
all nonzero component of Pi t j m (with i / j / m) are 
equal to j3. 

The solution of Eq. (|l^) under substitution of 
Eqs.( ^0t^2| ) can be written in the form 



2ttv 



i 



(14) 



where x = P \J / (ec^) is the dimensionless parameter 
of the piezoelectric interaction ( \ is the small parameter 
of our consideration), 7 g is the structural factor, which 
is determined by the boundary conditions on ip and by 
the type of the elastic mode. By using Eqs. fl[i|,[l4j), we 
arrive at the expression 



9qk 



2n , , 

—ev 1 Xlq(Zk), 



(15) 



where Zk is the z coordinate of the electron layer. The 
explicit expressions for r y q {zk) will be presented in Sec. IV, 
when we specify the geometry of the model. 



2m 



■(*V< 



1 



(18) 



where a M is the potential of the Chern-Simons gauge field, 
is the electromagnetic field potential, is fully asym- 
metric unit tensor, fip is the chemical potential, ip is the 
parameter of the model, which corresponds to the num- 
ber of the gauge field flux quanta carrying by the compos- 
ite quasiparticle (?/> is the even number), bi is an auxiliary 
field, which is introduced for modelling the Fermi liquid 
interaction. In Eq. ( |l8| ) the transverse gauge for the a 
field is implied. 

The value of £ can be determined from the condition of 
that the Lagrangian (|l^) provides the Galilean invariant 
value for the electric current (see, for example, Ref. |P4| ): 



SL 

JT' J '' 



en 



< b, >-- 



m b 



(19) 



where jt = — (e/m*)**(iVj + a L + eAj)*, n is the uni- 
form electron concentration. The mean value of the field 
b is determined by the equation 

< ^ >= = - < ji > < h > +£ < bi > . (20) 
obi e m* 



III. RESPONSE FUNCTIONS IN THE 
COMPOSITE FERMION MODEL 

Let us discuss the calculation of the response functions 
D for the composite fermion model. At the beginning, 
consider the monolayer electron system. The random 
phase approximation gives the answer 



D=(1-D V)- 1 D , 



(16) 



where V = 2ire 2 /eq is the Fourier component of the 
Coulomb interaction, and Do can be expressed though 
the polarization operator (ILjy) of the electromagnetic 
field 



Do — — tIIoo- 



(17) 



We use the approach which is the Lagrange formulation 
of the modified random phase approximation (Ref. 23 1) 
for the calculation of iloo- In Ref. [@ the Landau Fermi 
liquid interaction was taken into account. The F\ con- 
stant (see, for example, Ref 0)) °f that interaction is 
determined by the ratio m*/mb, where m* is the com- 
posite fermion effective mass, m b - is the band mass of 
electrons. It was shown in Ref. p3| ], that taking into ac- 
count that interaction one obtains the response function, 
which satisfy the /-sum rule. 

The Lagrangian of composite fermions can be written 
in the form 



The condition of consistency of Eqs. (|19|j20j) yields £ = 
n /(m* - m b ). 

The self-consistent effective field acting on the compos- 
ite fermions is B e ff — B(l — vip) (where B is the exter- 
nal magnetic field). At filling factors Vj — N/(Nip ± 1) 
(N is an integer value) B e ^ = ztnoipo/N, where (fo is 

the magnetic field flux quantum. The field B^. , cor- 
responds to the integer value N of filled Landau lev- 
els. Therefore, vfs correspond to the incompressible 
states. The field B** , defines the effective cyclotron en- 
ergy co ce = 2irno / m* N and the effective magnetic length 
A e// = {N/2nn )^ 2 . 

To determine the value of Iloo one can calculate the 
quadratic part of the effective action for the electromag- 
netic field. To do that we calculate the functional inte- 
gral over the fields \& in the expression for the partition 
function of the system, described by the Lagrangian (|l^) . 
Then we calculate the functional integral over the fields a 
and b in the vicinity of the saddle point. Here we present 
the final result at T = 



floo = 



? ? 
e q 



S 



2iru> ce Ai ' 



where 



So = £ -p(So(£ 2 + AO-£i), 
Ax = (1-VS 1 ) 2 -^E (E 2 + AT) 

F = Z 2 + N +(Lu/Lu ce ) 2 S 



pF, 



(21) 



(22) 
(23) 
(24) 



L = , i'*(idt + hf — «o — eAo — 



with p — (to* — m b )/{m*N), 
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N-l 



E i = (sgn(B eff )ye-* £ £ 



i™ ™ 1 (m — n) 
(uj/uj ce ) 2 - (m - n) : 



n=0 m=N 



jrm — nf \ 

x[(m-n-x)L™- n (x)+2x P- (25) 

In Eq. @ x = (qX eff ) 2 /2, L™~ n (x) is the generalized 
Laguerre polynomial. 

Substituting Eqs. (§J|-[r7]j2l]) into Eq. (|) we obtain 



&V 2 2/ \ fq S 



(26) 



where /, = e 2 q/(eiv ce ) 

The case of the double layer Hall system is considered 
by a similar way. Two types of the Chern-Simons fields 
and the fermion fields are introduced (they correspond 
to the two electron layers). The free part of the Chern- 
Simons Lagrangian includes non-diagonal terms with re- 
spect to the layer index. 



^cs — -^—tijaokMkk'diajk' 

47T 



with 



M = 



tp — s 
—s ip 



(27) 



(28) 



where s (an integer number) is the number of the Chern- 
Simons flux quanta of the type " 1" , carrying by the quasi- 
particle of the type "2" and vice versa (there exists the 
interlayer statistical interaction in the system). Here we 
specify the case of two equivalent electron layers. At 
ij) = s it is better to rewrite the Lagrangian in terms of 
the in-phase and out-of-phase combinations of the fields 
afe- 

The incompressible states for the double layer system 
correspond to the filling factors 



N 



N(ip + s) ± 1 



(29) 



(calculated for one layer). The expressions for B e ^, ui ce , 
\ e ff remain unchanged. Eq. ( |l6| ) obeys the matrix form 
with 



V = 



2we 2 
e<7 



1 e~ qc 
-qd i 



(30) 



where do is the distance between the electron layers. 

The expression for the velocity shift is modified to the 
following one 



_ 1 2 o , r (7 g (zi)+7 g (^)) 2 
v - 2 XboM A + -f q E + S 



+ 



(lq(zi) -lq{z2)Y 
A_ - f q E_S 



(31) 



where E± = 1 ± exp(— qdo) 



A ± = (l-(^±s)E 1 ) 2 

-(V-± s ) 2 So(s 2 + ^)-^. 



(32) 



The ground state of the double layer system is de- 
scribed by the quantum numbers i\) and s. One can see 
from Eq. (p^), that the same v may correspond to differ- 
ent sets of i/j and s (different phases). If the correlation 
between the layers is a small one (at large do), the phase 
with s — is realized. Reduction of the interlayer dis- 
tance may result in a phase transition into the state with 
s ^ (see, for instance, Ref. |2f|). The transition may 
be accompanied by the shift of the velocity of the acous- 
tic wave, interacting with the electrons, as described by 
Eq. @. 



IV. VELOCITY SHIFT FOR THE WAVEGUIDE 
ELASTIC MODES 



Eqs. (^6 31) define the value of Av relative the velocity 
for the systems without the electron-phonon interaction. 
In experiments the value of Av is measured relative to 
the velocity at B = 0. To compare the theoretical value 
with the experimental data one should know the value 
of (T X x{<l,u) (longitudinal component of the conductiv- 
ity tensor) at zero magnetic field. One can use the ap- 
proximation a xx — * oo at B = 0, but in that case the 
momentum dependence a xx (q,u>) at B — is neglected. 
We propose to measure the difference between the sound 
velocities at v = Vf and v = 1. The case of v = 1 is 
also described by Eqs. (^q), j3l|), if one sets ip = s = 0, 
N = 1 m* = mt>. 

Thus, the value we calculate is given by the expression 



Av Av 

v v 



Av 

v 



-Vf 



(33) 



Let us consider the structure AIq^Gclo^As — GaAs — 
AIq sGao/fAs, which incorporates two electron layers at 
Zk = ±et (do = 2a). The parameters (3 and e are implied 
to be uniform for the whole system. Then, the boundary 
conditions for ip. require the continuity of the values of 
tp(q, z) and d<p(q, z)/dz at the interfaces. 

The solution of Eq. @ for the even mode (Eq. (§)) 
gives an odd function of f(z), while for the odd mode 
(Eq. (Q)) - an even function of f(z). 

At v — 1/5 one can expect the phase transition be- 
tween (ip = 4, s = 0) and (ip = 2, s = 2) states. We see 
from Eq. (|3l|), that the transition is accompanied by the 
shift of the velocity of the even elastic mode, while the 
odd mode does not show such a behavior. 

For the even mode the values J q (zk) are found to be 



lq{a) = -lq(-a) = ■ 



_l cos(i]qa) 
1 + V 2 
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2(1 + A) 



-A 



(34) 



The dependencies Av/v on the wave vector are shown 
in Fig.l at -0 = 4, s — and ip = 2, s — 2. The param- 
eters (3 = 4.5 • 10 4 dyn^/cm, e=12.5, n = 10 11 cm" 2 , 
m& = 0.07m e (m e is the bare electron mass), to* = 4m;,, 
C44 = 6 • 10 11 dyn/cm 2 , p;=5.3 g/cm 3 , q,/c; = 1.05, 
do = 300A are used for the calculations. We see, that 
at finite q the value of the velocity shift is within the ex- 
perimental resolution, while it is rather small due to the 
smallness of Ac = (q, — c/)/c/. 




FIG. 1. Velocity shifts of the lowest even mode vs the wave 
vector for the double layer electron system at Cb/ci — 1.05. 1 
s = 2. q is in ^/4nno units. 



0; 2 



Now we analyze the geometry for which one could ex- 
pect larger values of the effect. Let the heterostructure 
is placed into the matrix of another material for which 
the value of Ac is much larger then for the previous case. 
For simplicity, we consider the heterostructure as a uni- 
form acoustic medium and assume that there are thin 
screening layers on the interfaces (then the dielectric and 
piezoelectric properties of the bulk matrix are unimpor- 
tant). The boundary conditions are tp(z = ±a) = 0. 

If the structure contains a single electron layer at z = 0, 
then the velocity shift can emerge only for the odd mode. 
The effect has a threshold character (at small q the odd 
mode does not exist) 

The structural factor in Eq. (Eq) is equal to 



7,(0) 



1 + 7f 



cosh(qa) / 



(35) 



The dependencies Av/v on q are shown in Fig. 2 at 
u f = 2/5,3/7,4/9 (N = 2,3,4) for the lowest odd elas- 
tic mode. The parameters Cb/q = 1.5, 2a = 150(M are 
used (other parameters are the same ones) . We see, that 
the function of Av(q) oscillates and the period of oscil- 
lation becomes smaller when Vf approaches to v = 1/2 
(it reflects the fact that the effective magnetic length is 
increased). 



Av/v, 
10" 4 




FIG. 2. Velocity shifts of the lowest odd mode vs the wave 
vector for the monolayer electron system at cj/c; = 1.5. 1 - 
v = 2/5; 2 - v = 3/7; 3 - v = 4/9; q is in ^/Anno units. 




FIG. 3. Velocity shifts of the third even mode vs the wave 
vector for the double layer electron system at Cf,/q = 1.5. 1 - 
s = 2. q is in \/47mo units. 



0; 2 



The absolute value of the velocity shift is much higher 
that for the previous case. The effect is lowered at larger 
a. 

If the heterostructure contains two electron layers at 
z = ±c?o/2, then the phase transition between ip = 4, 
s = and ijj = 2, s — 2 states is accompanied by the shift 
of the velocity of the even mode only, as for the previous 
case. The maximal value of the effect takes place for one 
of higher harmonics (we specify the case of do <C 2a). 

The structural factor is equal to 



7 9 (d /2) = -7,(-do/2) 

sm(jiqa) sinh(gdo/2) 
cosh(ga) 



n 



1 + r/ 
sin(77gd /2) 



(36) 



The dependencies of Av/v vs q are shown in Fig. 3 for 
the third even harmonics (do — 300A, 2a — 1500A). 
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Thus, the measurements of the velocity shifts for the 
non-uniform plane elastic modes of the waveguide type 
propagating in heterostructures AlGaAs — GaAs can be 
considered as a possible way of studying the dynamical 
properties of the monolayer and double layer quantum 
Hall systems. Unlike the case of the surface acoustic 
wave propagation, the restriction on maximal values of 
the wave vectors, for which the velocity shift can be ob- 
served experimentally, does not take place for the waveg- 
uide modes. In case of study of incompressible fractional 
Hall states, the effect can be used for measuring the 
dependence of the effective magnetic length on the fill- 
ing factor and for observing the phase transitions in the 
double-layer systems under the interlayer distance varia- 
tion 
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